O NMR study of the intrinsic magnetic susceptibility and spin dynamics of the 
quantum kagome antiferromagnet ZnCu 3 (OH) 6 Cl2 
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We report through 17 NMR, an unambiguous local determination of the intrinsic kagome lattice 
spin susceptibility as well as that created around non-magnetic defects arising from natural Zn/ Cu 
exchange in the 5=1/2 (Cu 2+ ) herbertsmithite ZnCu3(OH)eCl2 compound. The issue of a singlet - 
triplet gap is addressed. The magnetic response around a defect is found to markedly differ from 
that observed in non-frustrated antiferromagnets. Finally, we discuss our relaxation measurements 
in the light of Cu and CI NMR data and suggest a flat q-dependence of the excitations. 

PACS numbers: 



When nearest neighbor interacting 1/2 spins decorate 
the vertices of corner sharing triangles, as in the kagome 
antiferromagnet (AF), the geometry of the lattice effi- 
ciently frustrates the magnetic interactions and the clas- 
sical Neel states are certainly destabilized. Since Ander- 
son's seminal idea of a resonating valence bond spin liquid 
state [l[, the actual ground state has been a matter of 
deep conceptual debates and various theoretical propos- 
als were made in the last 15 years @, S S E, H, 0, @, S E3] 
with no definitive conclusion so far. The experimen- 
tal situation is quite contrasted and has been lagging 
well behind theory for long, unable to clarify the pend- 
ing debates on the nature of the kagome antiferromagnet 
ground state. The main reason is that achieving well 
decoupled perfect kagome planes has represented a real 
challenge since the pioneering work on kagome bilayers 
in SrCr 8 Ga 4 0i 9 [ll|. 

Quite recently the end member [x = 1) of the parat- 
acamite family Zn^Cu^n^OH^C^, known as herbert- 
smithite, under the name of the discovered mineral, was 
revealed as a "structurally perfect 5=1/2 kagome antifer- 
romagnet" [12} and opened a new avenue to comparisons 
between experiment and theory. Due to a more favor- 
able electrostatic environment, Cu 2+ is expected to pref- 
erentially occupy the distorted octahedral kagome sites. 
Therefore ZnCu3(OH) 6 Cl2 could feature perfect Cu 2+ , 
5=1/2, kagome planes separated by non magnetic Zn 2+ 
layers (Fig. 1). The lack of observation of any order down 
to the lowest explored temperature (50 mK~ J/3000 - J 
is the exchange interaction) [l3|, 14, HI is a clear indi- 



cation that this material is the best to-date to approach 
the ideal physics of a kagome lattice. 

Since the very first studies, a debate has settled about 
the actual low-T susceptibility x, one of the most crucial 
tests of the low-T ground state. Would a singlet ground 
state occur, \ should yield the singlet-triplet gap. In 
complete contradiction with this picture, the macroscopic 
susceptibility Xmacro, as measured by SQUID, shows a 



monotonous increase down to T =0 [lj| . Although addi- 
tional Dzyaloshinski-Moriya (DM) interactions were pro- 
posed to explain such a low-T increase [l6j], recent neu- 
tron structural refined data rather indicate 5-10% 
Zn 2+ /Cu 2+ inter-site exchange which could be responsi- 
ble for a Schottky-like field dependent anomaly detected 
in specific heat [18| and be the source of all or part of 
the low-T increase of Xmacro [ijj. This would result 
from both isolated Cu 2+ spins on the Zn 2+ site which are 
weakly coupled to the kagome layers and non-magnetic 
Zn 2+ defects in the Cu 2+ kagome planes which are known 
to generate a specific spin texture in their vicinity as was 
suggested in former studies of magnetic dilution of frus- 
trated lattices [HI, 0, EH . 

Local probes, such as used in nuclear magnetic reso- 
nance (NMR) , are invaluable to sort out the various con- 
tributions to X- They enable one to map out a histogram 
of susceptibilities at various sites, otherwise summed up 
in Xmacro- 17 NMR of enriched samples appears as the 
best compromise for such a local study. Due to the strong 
magnetic Cu-Cu bond through the bridging oxygen, one 
can indeed safely expect that O dominantly couples to 
the 2 neighboring Cu sites of the kagome planes only 
(Fig. 3) while its detection is not wiped-out through re- 
laxation effects at variance with the case of Cu [2J] . For 
the first time, we are able to track accurately the kagome 
lattice susceptibility from 300 to 0.47 K, and also, beyond 
establishing definitely the existence of Zn/Cu substitu- 
tion, use the latter perturbation to reveal some original 
facets of the kagome physics around these defects. We 
also present a comparison of our 17 NMR relaxation 
data with that from CI and Cu sites 24| , which suggests 



a non-dispersive branch of low energy excitations. 

The NMR experiments were performed on several un- 
aligned powder samples of different origins which gave 
identical results. Preparation and sample characteriza- 
tion follow standard routes described elsewhere [H, [l3| 
except that 50% 17 enriched water was used for the 
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preparation and samples were not washed after filtration. 
17 NMR spectra were recorded from room-T down to 
0.47 K by integrating the echo recorded by sweeping the 
field at a fixed frequency of 38.974 MHz. Delays between 
pulses were 10-20 fj,s, except in contrast experiments per- 
formed at low-T where delays, up to 200 /js allow to single 
out slowly relaxing sites. Quadrupolar and magnetic con- 
tributions were separated by taking additional spectra at 
19.974 and 29.974 MHz. Quadrupolar singularities typi- 
cal of a powder pattern are easily observed on the wings 
of our spectra at high-T (Fig. 1) and both a quadrupolar 
frequency vq = 0.70(3) MHz and the asymmetry param- 
eter r] = 0.64(10) were determined by the extreme cut- 
offs in the powder pattern ( 17 I = 5/2). The T\ recovery 
curves were obtained by irradiating the main line with 
a train of 5 pulses separated by ~ T%. Each component 
mixes all nuclear transitions (powder pattern) and we 
fixed the coefficients of the recovery law at the T = 70 K 
values. 

In the present letter, we track the evolution of the cen- 
tral part of the NMR spectrum. In the ideal case of a 
fully occupied Cu 2+ kagome plane, the structure yields 
a single environment of the O site therefore one should 
observe a unique central —1/2 — > 1/2 sharp singularity 
in the high-T range (T > J), the shift of which reflects 
the susceptibility of the two neighboring Cu. On the 
contrary, two different lines rather than one, correspond- 
ing to central transitions are surprisingly evident in our 
T = 300 K spectrum [(M) and (D) in Fig. 1]. This clearly 
indicates the existence of two distinct local magnetic en- 
vironments. The ratio of the magnetic shifts ~ 2 for the 
two lines (M) and (D) naturally supports the association 
of the most shifted and most intense line, with the vast 
majority of O nuclei locally coupled to two magnetic Cu, 
while for the other line (D), O couples to only one Cu, 
hence the other neighboring site is occupied by Zn. Since 
for each Zn at a Cu site of the kagome plane four O are 
affected, the 20(5)% intensity of this line is found consis- 
tent with a 6(2)% Zn/Cu substitution suggested by other 
techniques [3, El EE El ■ 

The shift of both lines, and K W can be easily 

tracked in a broad T-range. The main line (M) gives ac- 
cess to the intrinsic susceptibility of the kagome lattice 
through the line position while the physics of a Cu 1st 
near neighbor of a non-magnetic defect is best probed 
through the other line (D). The positions of the (M) and 
(D) lines were found to scale with the field after correc- 
tion of second order quadrupolar effects for low frequen- 
cies, which proves their magnetic origin. No simulation 
even allowing different axis for the EFG and shift ten- 
sors could account for the presence of these two peaks in 
a single O site model. 

The (M) line is found to shift towards low fields (in- 
crease of susceptibility) between 300 and 100 K (Fig. 2a), 
goes through a maximum between 100 and 50 K and 
then progressively moves back towards the reference be- 
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FIG. 1: 17 NMR spectra from 300 K to 0.47 K. The shift 
reference is indicated by the vertical line. The line on the 
175 K spectrum is a simulation using two components, (M, 
dotted line) and (D, hatched), see text. Large full circles 
indicate the position of the center of the line which would be 
expected from Xmacro- 

low 40 K while a rapid increase of the width occurs 
(Fig. 1). In the low-T regime, below 5 K, the linewidth 
tends to saturate and no difference was found between 
spectra taken at 1.3 and 0.47 K. The variation of the 
shift, K^ M \ is reported in Fig. 2b. It monitors the 
T-evolution of the susceptibility of Cu in the kagome 
planes, except those sitting nearby non-magnetic defects, 
through K^ M ' = 2 Ay On , with the hyperfine constant 
A = 35(2) kOe//iB [25j |. It is worth noting that this cou- 
pling value is approximately 30 times larger than those 
reported for CI NMR and /iSR experiments, hence the 
invaluable sensitivity of the 17 O probe. 

Two main model-free results can be deduced from 
K( M ); (i) A maximum of the kagome susceptibility is 
found to occur around J/2 which differs markedly from 
Xmacro- Defects induced by Cu/Zn substitution and re- 
sponsible for the low T increase of Xmacro, contribute to 
the NMR linewidth. (ii) The susceptibility reaches a non- 
zero value at low T, typically 1/3 of the maximum value 
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FIG. 2: Top: blow-up of high T spectra. The shift of the main 
and defect lines are evident and found different. Bottom: Plot 
of the magnetic shift of both lines versus T. The line through 
the K (D) plot represents K {M) /2. The sketch on lower left 
corner is for what would be expected around a Zn substituted 
on the Cu kagome plane at T = and thick lines represent 
Cu-Cu dimers. 



which indicates the absence of a gap in herbertsmithite. 

One can invoke three possible scenarii: (i) a non-singlet 
ground state of the ideal kagome system, such as observed 
here, with a decreasing susceptibility at low T associ- 
ated with the strengthening of the AF correlations below 
T ~ J/2 which then level off as T — > 0, and surpris- 
ingly, do not yield a gapped ground state, (ii) Sticking 
to a singlet-triplet gap: The finite T = susceptibility 
observed in herbertsmithite may result from a filling of 
the hypothetical gap -expected to be at most J/20- un- 
der the modest applied field 6.5 T of our experiment. It 
is however not clear to us whether the magnetic energy 
brought to the system ~ J/30 should be compared to 
the gap value or to the larger one at which a magnetiza- 
tion plateau should occur. In addition, the value of the 
maximum of K^ M ^ around 50 K is not on-line with the ex- 
pected small gap value, (iii) The existence of additional 
DM interaction H(^x S^) between Cu 2 + spins which 
results from the absence of inversion center between Cu, 
as earlier suggested This would modify the T=0 

susceptibility by mixing singlet and triplet states. Al- 



though x could not be calculated below 0.25 J |l6[, we 
anticipate from a comparison of our K^ M ^ curve and fill ] 
that this DM term might be smaller than suggested and 
rather less than 0.05 J at least for one component of D. 
Finally, one should point that exact calculations, reliable 
down to 0.2-0.3 J, do not indicate any maximum of x- 
The deviance of our experimental data from the ideal 
case could be attributed either to the DM interaction or 
to the rather strong dilution of the network favoring a 
dimerization of the ground state (see below). 

We now turn to the defect line (D) which reflects the 
susceptibility of a Cu in the immediate vicinity of a non- 
magnetic Zn substituted on a kagome site. At high-T, the 
(D) line follows a trend similar to the (M) line, as can 
be observed through the comparison of spectra taken be- 
tween 300 and 150 K (Fig. 2a). The 1:2 scaling of 
with R( M ' (Fig. 2b) indicates that in this T-range defects 
have negligible effects on the Cu sites in their immediate 
vicinity, which underlines the weakness of the correlation 
length even for T ~ J, another well-known landmark 
of frustration. The (D) line surprisingly appears below 
50 K as a sharp prominent feature of the overall spec- 
trum, which contrasts with the broadening of the (M) line 
in the same T-range (Fig. 1). This relative narrowness 
points at a reduced distribution of susceptibilities which 
witnesses a quite unique environment around the defect. 
In addition, below 50 K, the decrease of the shift is found 
sharper than for the (M) line, an indication that the sus- 
ceptibility of Cu in the neighborhood of a Zn decreases 
faster. Finally, the transverse relaxation was found sub- 
stantially longer than for the (M) line in this T-range 
which enabled us to perform contrast experiments and 
better isolate the defect line. A relative intensity of the 
(D) line ~ 30(4)% can then be calculated. 

The picture which can be sketched differs markedly 
from the situation met in non-frustrated AF where the 
spin response is commonly enhanced in the immediate 
vicinity of a non-magnetic defect [26| and would yield 
a rapid low-T increase of K^ D ' rather than a decrease 
as clearly observed here. Our findings arc in qualitative 
agreement with the results from exact diagonalizations 
where a defect is found to induce a quite reduced suscep- 
tibility or even singlets on the 4 neighboring sites while an 
enhanced stag gered susceptibility is induced for further 
neighbors [22j . The latter would be naturally responsi- 
ble for the smoothing out of the quadrupolar singulari- 
ties and could explain part or all of the broadening of the 
(M) line. Another source of broadening could also result 
from the alternating DM contribution as proposed re- 
cently [l6| • In the ideal case of singlets induced around a 
non-magnetic defect, a simple counting argument yields 6 
O sites per defect (fig. 3) with a zero shift, therefore 36% 
of the O sites, in agreement with the detected intensity. 
The decrease of K^ D \ sharper than K^ M \ could there- 
fore be explained by a tendency for Cu to dimerize near 
a defect, as expected from the local relief of frustration. 
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FIG. 3: Semi-logarithmic T 1 1 plots for 17 (this work) and 
for 63 Cu and 35 CI HI). The similarity of the T- variation for 
the 3 nuclei is evident. Inset: 17 T 1 _1 versus T -1 . A non-zero 
gap would yield a linear variation, not observed here. 



imum in the susceptibility at ~ J/2 certainly related to 
the kagome physics and also observed in other corner- 
sharing kagome based lattices [3, Hi], 27 1 . A rapid de- 
crease of the susceptibility to a finite value at low-T 
and ungapped excitations point at a non-singlet ground 
state in Herbertsmithite. The excitations have not a 
marked structure in g-space, a landmark of frustration. 
A depressed susceptibility is observed for Cu next to de- 
fects. This response to defects is certainly specific to the 
kagome physics and observed for the first time here. It 
opens the route for revealing intrinsic physics through 
magnetic perturbations after a better control of Zn/Cu 
substitution is achieved. In order to clarify whether a 
gap would occur in ideal kagome Heisenberg antiferro- 
magnets, an accurate experimental determination of the 
DM term is necessary and the search for a zero-defect 
" second generation" compound guides future experimen- 
tal work. 



The fact that we do not observe a zero-susceptibility and 
a gap resulting from dimers could certainly be ascribed 
again to DM interactions. 

Finally, T\ measurements performed at the (M) peak 
yield information about the excitations. A second con- 
tribution, likely due to the proximity of the (D) line was 
found below 50 K, with a T\ which gets longer when T 
is lowered. This certainly corroborates the decrease of 
magnetic effects in the vicinity of a Zn defect, as already 
suggested by our shift data. In Fig. 4, we plot the short 
contribution below 70 K, dominated by the physics of 
the kagome plane, and compare it to existing CI and O 
data 24j. Between 70 and 30 K, we do not find a marked 
maximum such as observed for CI NMR which had been 
suggested to be of structural origin. Due to its weak 
quadrupolar moment and its larger coupling constant to 
Cu, 17 NMR is less sensitive to structural effects and 
probes dominantly magnetic properties. Whatever the 
origin of the 35 (1/Ti) peak, our data indicates that mag- 
netic properties are not affected. Like for Cu T^ 1 , we 
find a downward curvature of our data below 20 K and 
the inset of Fig. 3 clearly shows that the excitations are 
not gapped. The absence of anomaly in the T\ variation 
is on-line with the absence of signature of a transition 
as observed in our shift measurements down to 0.47 K 
or in earlier measurements 
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|. We rather find 
that Tj - ~ T°- 73 ( 5 ) and the perfect agreement of all data 
taken on various probes with different local g-space fil- 
tering of fluctuations is a clear experimental indication 
that the excitation modes responsible for this relaxation 
are dispersionless as suggested for long [29j]. One would 
indeed expect that "triangular 120°" modes or (it ,ir) 
modes are filtered out and do not contribute to the re- 
laxation respectively at the CI or O sites at variance with 
Cu. 

To summarize, we demonstrate the existence of a max- 
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